Momentum space anisotropy in doped Mott insulators 
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We study the single hole tt't" J -model numerically to address the momentum space anisotropy 
found in doped Mott insulators. A simple two band picture to understand how the doped hole is 
screened by the spin background in states of different momenta is proposed. In this picture, the 
disparity between the nodal and antinodal regions, observed by experiments in the underdoped 
cuprate superconductors, follows from local energetic considerations and amounts to the distinction 
between well defined quasiparticle behavior and spin-charge separation related phenomenology. 



Introduction - It is known, both experimentally and 
theoretically, that doped Mott insulators can be signifi- 
cantly anisotropic in momentum space. Explaining this 
behavior is relevant to understanding the properties of 
metals close to the Mott insulating state, like the pseu- 
dogap regime in high-Tc superconductors. P, Q 

The pseudogap measured by ARPES is the energy 
difference hetwe&a. the one-electron spectral features at 
the nodal [k = (±-|,±|-)] and antinodal points [k ~ 
(tt, 0), (0, tt)]. "i) This difference, which increases as both 
the temperature and the level of doping in the Cu02 
layers are reduced, is accompanied by the anisotropic 
behavior in fc-spacepresent in various experimental ob- 
servations. 0, HIllllEIHllil Indeed, the deeply 
hole underdoped regime shows a clear nodal-antinodal 
dichotomy as low energy quasiparticle peaks exist around 
the nodes while there is no evidence of quasiparticle-like 
behavior near the antinodal points. ^ S 13 In elec- 
tron underdoped cuprates, low energy spectral weight 
appears around (tt, 0) instead as the pseudogap is pushed 
toward the zone diagonal. ^10] Interestingly, Raman spec- 
troscopy and the violation of the Wiedemann-Franz 
law 12| in electron doped materials suggests the presence 
of chargeless excitations around (^, ^). Hence, the pseu- 
dogap in the one-particle dispersion appears to strongly 
reduce the electronic character of excitations in the pseu- 
dogap region. 

Strongly correlated materials close to the Mott insu- 
lator transition are notorious for the near degeneracy of 
different ordered states, some of which have been argued 
to provide a scattering mechanism that reduces quasipar- 
ticle features in the pseudogap region, 
this letter, we propose a different scenario to understand 
the observed momentum space anisotropy. In particular, 
we argue that states close to (f , f ) and (tt, 0) are dif- 
ferent because, due to the competition between the spin 
exchange energy and the hole kinetic energy, spins sur- 
rounding the hole in states near (|-, ^) behave differently 
from the spins surrounding the hole in antinodal states. 
(iM (31 The nodal-antinodal dichotomy then reflects the 
existence of two distinct ways in which the lattice spins 
screen the doped holes. 

The 2D tt'i"J-model, which is the simplest model 
to study doped Mott insulators as of relevance to the 



cuprate systems, reproduces the experimental pseudo- 
gap dispersion and accounts for anisotropic behavior in 
/c-space. ^15, 16, 17, 18, 19] Below, we employ the ex- 
act diagonalization (ED) and the self-consistent Born ap- 
proximation (SCBA) techniques to study the single hole 
ti'i" J-model and to address how the hole is screened by 
the local spins in different regions of momentum space. 
We find that the single hole states can be understood 
as the superposition of two distinct states, namely a 
hole-like quasiparticle state and a state where the hole 
strongly distorts the surrounding spins. Due to their dif- 
ferent properties, these states predominate in different 
parts of /c-space in the pseudogap regime, thus leading to 
the disparity between the nodal and antinodal regions. 
Our results are valid for both hole and electron doped 
materials - for our purposes, the main distinction be- 
tween the two regimes is that in the former the pseudo- 
gap opens in the antinodal region while in the latter it 
opens in the nodal region. 

Model - The single hole 2D tt't" J Hamiltonian is 
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equals t, t' and t" for first, second and third 



where t. 

nearest neighbor sites respectively and vanishes other- 
wise. Ci^cr is the constrained operator Ci^^ = Ci,cr(l — 
rii^^a-)- The exchange interaction only involves nearest 
neighbor spins for which — J. We only consider 
0.2 < J < 0.8 (units are set so that t — 1). Unless 
otherwise stated, all our results come from the ED of 
Htt't" J on a 4 X 4 lattice. Since we want to analyze how 
the hole affects the local configuration of the surround- 
ing spins we believe that the study of such a small lattice 
is relevant. We also present results from the SCBA ap- 
proach to the tJ model [13, mil on a 16 X 16 lattice 
to further support the ED analysis. 

We start by considering the t' ,t" = case. In partic- 
ular, we use the ED technique to determine the lowest 
energy state for each momentum, denoted by jf/jfe, J), for 
both J = 0.2 and J — 0.6. Quite surprisingly, we find 
that if we perform the same calculation for different val- 
ues of J the resulting states \tpk,J) have almost complete 
overlap with the Hilbert space {li/'fc,^/ = 0.2),|?/;fc, J — 
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TABLE I: Square of the overlap of \^k,J,t',t") with the 
Hilbert space {\Tpk,J = 0.2,t' ,t"),\i!k, J = 0.6,t',t")} for 
different J and fe. Both ED and SCBA results are shown for 
t',t" = 0. ED results are also shown for t' = -0.2, t" = 0.1. 
For t' = -0.2, t" = 0.1 and fe = (7r,0) the Hilbert space 
{\'tpk,J = 0.4,t',t"), \tpk,J = 0.8, t',t")} was used instead. 



0.6)} (Tabic HI) . The SCBA technique leads to the same 
conclusion (Table QJ, thus supporting that this result is 
not specific to the 4x4 lattice used in the ED calculation. 

Therefore, for 0.2 < J < 0.8, the wave function IV'fc, ^) 
can be written approximately as \ipk,J) — <l{J)\Qk) + 
u{J)\Uk)- Here, q{J) and u{J) are J-dependent co- 
efficients obeying the normalization condition q(JY + 
u{JY — 1- \Qk) and \Uk) are orthonormal states that 
span the 2D Hilbert space {|^/'fe,J = O.2),|'0fc,J = 
0.6)}. Of course, there is freedom in choosing \Qk) 
and \Uk)- A physically sensible choice, though, comes 
from requiring that g(J) increases with J while u(J) 
decreases. Since increasing J enhances the quasiparti- 
cle features of doped carriers, [i^ this is automatically 
satisfied if the quasiparticle spectral weight vanishes for 
\Uk), i-e. {Uk\ck,a\ilF GS) = where |HF GS) is the 
groundstate of the half filled system. This prescrip- 
tion uniquely determines Q states {\Qk)) and U states 
i\Uk)), |2^ which are not eigenstates of Hu't"j, as long 
as (V-fc, J|?fc,.|HF GS) ^0. 

For general k and t',t" the state \^k, J,t' ,t") only 
has non-zero quasiparticle spectral weight for J larger 
than a certain Jc{k,t',t"). The above prescription 
can be applied for J > Jc{k,t',t"), in which case 
it is a good approximation to write {ipk, J,t' ,t") = 
q{J,t\t")\Qk,t',t") + u{J,t',t")\Uk,t',t"). This fact is 
illustrated in Table U for t' = -0.2, t" = 0.1. Note that 
0.3 < Jc < 0.4 for k = (7r,0), in which case we used 
the eigenstates for J = 0.4 and J = 0.8 to determine 
IQfe=(7r,o): *") and \Uk=(7,.o),t' ,t")- As discussed be- 
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FIG. 1; Average value of the staggered magnetization on 
sites that are z/"* nearest neighbors to the hole, (i^) = 

and S^{,.) = ((-)'-+'" 5fe^(i))., for differ- 
ent momenta fe. (a) and (b) t' — — 0.3, t" — 0.2. (c) and (d) 
t',t" = 0. (e) and (f) t' = 0.3, t" = -0.2. 



low, for J < Jc{k,t' ,t") the eigenstate \tpk, J,t' ,t") dis- 
plays properties akin to those of U states and we denote 
it by \Uk, J, t' , t") (where the tilde is used to distinguish 
it from the U state derived for J > Jc{k,t' ,t"), which 
has no J dependence). 

The previous argument points out that the eigenstates 
{ipk, J,t't," ) have features that get enhanced by an in- 
creasing J and features that become more pronounced 
when J is reduced. By definition, Q and U states capture 
these features and, not surprisingly, below we show they 
display distinct physical properties: in Q states the spin 
antiferromagnetic (AF) order is robust to the presence of 
the hole while in U states the spins rearrange to facilitate 
hole hopping. On general grounds, a construction similar 
to the above one might be valid for models other than the 
2D tt't" J-model. The significant fact about this model is 
that, for experimentally relevant parameters, the overlap 
of both Q and U states with \il!k,J,t', t") is large and ex- 
hibits a considerable momentum dependence which leads 
to the anisotropy in /c-space. ji^l 

Properties of Q and U states - To explore how the 
hole affects the spin background in Q and U states con- 
sider the spin density pattern around the hole, (i) = 

(>^=IE,'5|+,s-,+i/2c]_+i/2l^fe) + N^A With Y = Q,U, 
[i^j l as well as the hole momentum distribution function. 



'{q,a) = {Yk\c. 



\Yk) for Y = Q,U. Figs. 



Ma),^c) and^^e) depict the staggered magnetization in 
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TABLE II: E;4(q,-5) = E;(t^fc|' 



1 c 



l\Uk). q = 



(0,0) results involve sum over q = (0,0), q = (±|-,0) and 
q = (0, ±-|). q = (7i",7r) results involve sum over q = (tt, tt), 
q = (±f ,7r) and q = (7r,±f ). J = 0.4. 



Q states, (-)'^+''' S'^(i), for different fe, i' and t". Since 
these states have a well defined quasiparticle character, 
the doped hole coexists with the AF spin pattern inher- 
ited from the undoped system. For the same reason, 
the hole momentum distribution function n^{Q,^^) is 
peaked at q = k while a smaller peak is also observed at 
q = fe + (tt, tt) due to the strong AF correlations. [2^ 

U states have no quasiparticle weight and display quite 
different behavior. Indeed, Figs. ^b),njd) andQJf) show 
that in \Uk) the AF spin pattern of the undoped system 
is destroyed and the staggered magnetization around the 
hole [given by (— '5'^(i)] is very close to zero and 
even negative. Moreover, the hole momentum distribu- 
tion function 7i^(q, — i) peaks around q — (tt, tt) for all 
momenta fc [Fig. |21 (a)]. The same real and momen- 
tum space properties were checked to hold for U states 
[these are the energy eigenf unctions J,t' ,t") when 
J < Jc{k,t' ,t")]. This fact is illustrated for U states 
concerning the parameter regime relevant to both hole 
doped cuprates (J = 0.4, t' — —0.3, t" — 0.2) and elec- 
tron doped cuprates (J = 0.4, t' = 0.3, t" = -0.2) [l3 
in Table In] which explicitly shows that in these states 
the hole density also peaks around (tt, tt) independently 
of the momenta fc. 

According to the above spin density results the ex- 
tra S"^ = — i spin introduced by doping spreads away 
from the vacancy in both U and U states. The resulting 
loss of spin exchange energy is accompanied by a gain in 
hole kinetic energy. Indeed, the hole momentum distribu- 
tion results support that, in these states, the hole always 
lies around the bare band bottom [which is located at 
(tTjTt)]. This evidence resembles predictions from spin- 
charge separation scenarios. Indeed, within the slave- 
boson formalism, the electron decays into a charged spin- 
less holon, which condenses at (tt, tt), and a chargeless 
spinon, which describes the delocalized spin-^ that car- 
ries the remaining momentum. We should remark that 
our calculation involves equal time properties in a small 
lattice and does not aim to prove the existence of true 
spin-charge separation. However, it supports that in U 
and U states the lattice spins screen the hole in confor- 
mity with short range aspects of spin-charge separation 
phenomenology. 
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FIG. 2: (a) Eq 'Jfe (q, ^|)- Empty symbols involve sum over 
q — (7r,7r), q — (±-|,7r) and q — (7r,±-|). Full symbols 
involve sum over q = (0,0), q = (±f ,0) and q = (0, ±-|). 
(b)-(d) Dispersion relations for \Qk) (full line), \Uk} (dashed 
line) and \tpk} (dash-dot line). Upper, middle and lower set 
of dispersions are obtained for J equal to 0.2, 0.4 and 0.7 
respectively, (o) indicates the best energy obtained by a linear 
combination of \Qk) and \Uk) when J < Jc{k,t' ,t") (in which 
case \tPk) = \Uk,J,t',t")). 
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TABLE III: AE'^, AE^ , AE"* and with k = k' = (tt, 0) 
and fc = fc" = (f , f ) for several t' and t" and J = 0.4. 



The previous results demonstrate that the local con- 
figuration of spins encircling the hole is quite different 
in Q and U states. In this context, it is important to 
remark that the influence of t',t" on the hole dispersion 
is sensitive to the surrounding spin environment. For in- 
stance, intrasublattice hopping is not frustrated by AF 
correlations and, indeed, the hole in Q states (which is 
surrounded by a spin configuration reminiscent of the un- 
doped AF groundstate) strongly disperses along the AF 
Brillouin zone boundary (AFBZB) for t',t" ^ (see Fig. 

The opposite limit occurs, for example, in certain 
U{1) spin liquids, where the underlying spin correlations 
inhibit coherent intrasublattice hopping. This limit 



4 



is closer to what is observed in U states, as t' ,t" are 
heavily renormalized by the spin background (whose AF 
correlations are depleted by the hole nearest-neighbor 
hopping processes). In fact, Table UTTl explicitly shows 
e9^ _ is almost one order of 



that A£:'5 = e9 - ^1 ,„ ,„ 

(tt.O) {7r/2,7r/2 

magnitude larger than Ai?^ 



-^(^/2,^/2)' where 



= {Qk\Hu't"j\Qk) and E'^ = {Uk\Hu't" jpk) ■ 
Momentum space anisotropy - In the tt't"J- 
model the intrasublatticc hopping parameters t',t" con- 
trol the dispersion i?^ = {ipk\Htt't"j\'ipk) along the AF- 
BZB and, thus, the pseudogap energy as well (this is the 
energy difference between the antinodal and nodal states 



AE'i' = e: 



(77,0) 



■E, 



■0 

(77/2,77/2) 



). Indeed, both Fig. Eland Ta- 



ble ^O] show that these parameters set the magnitude of 
the pseudogap, as well as its location in fc-space. [ifllE^ 
Interestingly t',t" also control the difference in the 
quasiparticle character of (7r,0) and (§, f ) states, as 
supported by the dependence of the overlap integral 
= \{^k\Qk)\'^ on these parameters fTable lllUl . The 
resulting t' , t" driven momentum space anisotropy is par- 
ticularly evident when both Q and U states have large 
overlaps with the eigenstates iV'fc), as it is the case for 
the value of J which is relevant to the cuprate systems, 
namely J = 0.4. [Ulii 

The main message of this letter is that the above roles 
of t' , t" can be understood with a simple two band picture 
involving Q and U states which captures the screening 
properties of the local spins encircling the hole. Indeed, 
as previously discussed, intrasublatticc hopping is frus- 
trated in U states and the energies E^^ and E^^^^ ^^2) 
have a small dependence on t',t". In contrast, t' < 
and t" > clearly increase the energy E'^^ and clearly 

decrease -E'(^/2 77/2) (^^&- El- Since It/jk) is a superpo- 
sition of \Qk) and \Uk), these t',t" induced changes in 
the energy of Q states push E^^ upward and lower 

-^(77/277/2)- Therefore, they increase the pseudogap en- 
ergy AE"^ . These changes also increase the energy differ- 
ence Sgo) 



■^(^,0) wl^ile reducing E'^^^ 



and, as a result, W^^. and 1^(^/2 7r/2) 
and larger respectively. For J 



^(77/2,77/2) 



"(77/2,77/2) 

become smaller 
0.4,t' = -0.3,t" = 0.2 
the energy E'^^ is so large that the minimum energy 
obtained by a linear combination of \Q{-„fi)) and |C/(77,o)) 
becomes higher than that of a different state |C/(77,o)). As 



a result, W, 



Q 



0. Since W, 



Q 



'(77,0) - "^'"'^^ "(77/2,77/2) - 0.75, the (|, |) 
state has strong quasiparticle character and a sharp dis- 
parity between nodal and antinodal states is encountered 
for these parameter values. Similar behavior is obtained 
when i' > and t" < with the role of momenta (tt, 0) 
and (|, f ) interchanged f Table llTT)l . 

The above argument shows that, in the pseudogap re- 
gion, the kinetic energy of a hole surrounded by an AF 
configuration of spins increases as the pseudogap gets 
larger. When this kinetic energy exceeds a certain value 



the spins around the hole lose their staggered pattern 
and, at short distances, the doped spin and charge sepa- 
rate. The resulting nodal-antinodal dichotomy, found to 
occur in the parameter regimes concerning hole and elec- 
tron doped cuprates, follows from local energetic consid- 
erations and not from the interaction of holes with other 
holes, an incipient order parameter or disorder. It also 
agrees with previous numerical evidence for spin-charge 
separation phenomenology at high energy, [la, llal 

Following our results, a new mean field theory of the 
tt'i" J-model was recently proposed that describes holes 
dressed by two different spin configurations. 30] It ac- 
counts for the observed momentum space anisotropy and 
is in good agreement with the doping dependence of 
ARPES experiments. 
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